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CN . We discuss how leptogenesis can explain the observed baryon asymmetry and sum- 

bJOj marize attempts of testing leptogenesis. We first perform estimates and discuss the 

main physics, and later outline the techniques that allow to perform precise compu- 
. tations. 
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1 Introduction 



The universe contains various relict particles: 7, e, p, v, He, Deuterium,. . . , plus likely some Dark 
Matter (DM). Their abundances are mostly understood, with the following main exception: 

nB ~ n B = M = (6.15 ± 0.25) • 10- 10 (1) 

where n 7 and n# and are the present number densities of photons and baryons (anti-baryons have 
negligible density). This is the problem of baryogenesis. Before addressing it, let us briefly summarize 
the analogous understood issues. 

As suggested by inflation, the total energy density equals the critical energy density, discussed later. 
Next, almost all relative abundances can be understood assuming that these particles are thermal relics 
of a hot Big-Bang phase. The number densities of electrons and protons are equal, n e = n p , because 
nothing violated electric charge. The relative proton/neutron abundancy was fixed by electroweak 
processes such as nv e <-> pe at T ~ few MeV. Neutrons get bound in nuclei at T ~ 0.1 MeV: the 
measured nuclear primordial nuclear abundances agree with predictions: n4 He /^p ~ 0.25/4, n-£,/n p ph 
3 10 _5 /2, etc. The neutrino density, predicted to be n Ve = n>v e = 3n 7 /22, is too low to be 
experimentally tested: the baryon asymmetry problem is more pressing than the analogous lepton 
asymmetry problem because we do not know how to measure the lepton asymmetry. Finally, the 
DM abundancy suggested by present data is obtained if DM particles are weakly-interacting thermal 
relicts with mass 

m ~ V^now • Mpi ~ TeV (2) 

where T now ~ 3 K is the present temperature and Mp\ ~ 1.2 10 19 GeV is the Planck mass. The 
LHC collider might soon produce DM particles and test this speculation. It is useful to digress and 
understand eq. (2), because the necessary tools will reappear, in a more complicated context, when 
discussing leptogenesis. 



1.1 The DM abundancy 

First, we need to compute the expansion rate H (t) of the universe as function of its energy density p{t). 
Let us study how a homogeneous p(t) made of non-relativistic matter evolves according to gravity. A 
test particle at distance R from us feels the Newton acceleration 

~ GM(R) ATrGpjt) 

R = = — R (3) 

where M (R) is the total mass inside a sphere of radius R and G = 1/Mpj is the Newton constant. By 
multiplying both sides of eq. (3) times R and integrating taking into account that p(t) oc l/R? (t) one 
obtains as usual the 'total energy' constant of motion, here named k: 



d 
di 



2 3 H 



TT o R 2 8vrG k 
sothat H 2 = — = -—p-— (4) 



Let us discuss the special case k = 0. It is obtained when the density p equals the 'critical density' p = 
Per = 3H 2 /8irG. k = is special because means zero 'total energy' (the negative gravitational potential 
energy compensates the positive matter energy): a universe with critical density that expands getting 
big for free could have been theoretically anticipated since 1687. Today we abandoned prejudices 
for a static universe, and more advanced theories put the above diseussion on firmer grounds. In 
general relativity eq. (4) holds for more general sources of energy (relativistic particles, cosmological 
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Figure 1: (a) Main reactions that determine primordial nuclear abundances. (b) How CMB 
anisotropies depend on the baryon abundancy Qb = Pb/Pct, compared with data. 

constant,...), and the constant k gets the physical meaning of curvature of the universe. The inflation 
mechanism generates a smooth universe with negligibly small k. 

Second, we need to know that a gas of particles in thermal equilibrium at temperature T>m has 
number density n eq ~ T 3 and energy density p cq ~ T 4 : one particle with energy ~ T per de-Broglie 
wavelength ~ l/T. The number density of non relativistic particles (T <C m) is suppressed by a 
Boltzmann factor, n eq ~ (mT) 3 / 2 e _m//T , and their energy density is p cq ~ mn DCl . 

We can now understand eq. (2), by studying what happens to a DM particle of mass m when the 
temperature T cools below m. 

Annihilations with cross section cr(DM DM — ► SM particles) 
try to maintain thermal equilibrium, tidm oc exp(— m/T). 
But they fail at T < m, when hdm is so small that the collision 
rate T experienced by a DM particle becomes smaller than the 
expansion rate H: 

T ~ n D Mff ^ H ~ T 2 /M P1 . 

As illustrated in the picture, annihilations become ineffective, 
leaving the following out-of-equilibrium relic abundancy of 
DM particles: 

»dm m 2 /M Pl o- 1 Pdm(T) 

; rs^ l.e. 

n 7 m 3 Mp\am P-y{T) 

Inserting the observed DM density, pdm ~ p 7 at T ~ T now , and a typical cross section a ~ g 4 /m 2 
gives eq. (2) for a DM particle with weak coupling g ~ 1. A precise computation can be done solving 
Boltzmann equations for DM. 

1.2 The baryon asymmetry 

Let us summarize how the value of the baryon asymmetry in eq. (1) is measured. The photon density 
directly follows from the measurement of the CMB temperature and from Bose-Einstein statistics: 
n 7 ~ T 3 . Counting baryons is more dimcult. Direct measurements are not accurate, because only 




m npM 
T n^y 



1 

M P iaT' 



(5) 
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some fraction of baryon formed stars and other luminous objects. Two different indirect probes point 
to the same baryon density, making the result trustable. Each one of the two probes would require a 
dedicated lesson: 

1. Big-Bang-Nucleosynthesis (BBN) predictions depend on ns/n 7 . Fig. la illustrates the main 
reactions. The important point is that the presence of many more photons than baryons delays 
BBN, mainly by enhancing the reaction pn <-» D7 in the <— direction, so that D forms not 
when the temperature equals the Deuterium binding energy B ~ 2 MeV, but later at T ~ 
B / \n{nB /n-y) ~ 0.1 MeV, giving more time to free neutrons to decay. A precise computation 
can be done solving Boltzmann equations for neutrino decoupling and nucleosynthesis. 

2. Measurements of CMB anisotropies [1], among many other things, allow us to probe acous- 
tic oscillations of the baryon/photon fluid happened around photon last scattering. A precise 
computation can be done evolving Boltzmann equations for anisotropies, assuming that they 
are generated by quantum fluctuations during inflation: fig. lb illustrates how the amount of 
anisotropies with angular scale ~ l/l depends on ub/uj. Acoustic oscillations have been seen 
also in matter inhomogeneities [1], at ~ 3<r level. 

2 Baryogenesis 

The small baryon asymmetry ns/n^ <C 1 can be obtained from a hot big-bang as the result of a small 
excess of baryons over anti-baryons. We would like to understand why, when at T < m p matter almost 
completely annihilated with anti- matter, we survived thanks to the 'almost': 

n B -n B (x 1000000001 - 1000000000 = 1. 

This might be the initial condition at the beginning of the big-bang, but it would be a surprisingly 
small excess. In inflationary models it is regarded as a surprisingly large excess, since inflation erases 
initial conditions. 

In absence of a baryon asymmetry an equal number of relic baryons and of anti-baryons survive 
to annihilations at T<m p . This process is analogous to DM annihilations studied in the previous 
section, so we can estimate the relic baryon density by inserting m = m p and a typical pp cross section 
a ~ 1/rrip in eq. (5), obtaining n p /n^ ~ m p /Mpi ~ 10~ 19 . Therefore this is a negligible contribution. 

Assuming that the hot-big-bang started with zero baryon asymmetry at some temperature T S> 
m p , can the baryon asymmetry can be generated dynamically in the subsequent evolution? Once that 
one realizes that this is an interesting issue (this was done by Sakharov), the answer is almost obvious: 
yes, provided that at some stage [3] 

1. baryon number B is violated; 

2. C and CP are violated (otherwise baryons and antibaryons behave in the same way); 

3. the universe was not in thermal equilibrium (we believe that CPT is conserved, so that particles 
and antiparticles have the same mass, and therefore in thermal equilibrium have the same 
abundance). 

Having discussed in section 1.1 a concrete out-of-equilibrium situation, it should be clear what the 
general concept means. 
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2.1 Baryogenesis in the SM? 



A large amount of theoretical and experimental work showed that, within the SM, Sakharov conditions 
are not fulfilled. At first sight one might guess that the only problem is 1.; in reality 2. and 3. are 
problematic. 

1. Within the SM B is violated in a non trivial way [3], thanks to quantum anomalies combined 
with extended SU(2)x, field-configurations: the anomalous B and L symmetry are violated, while 
B — L is a conserved anomaly-free symmetry. The basic equation is d^Jg ~ N gcn F® u F® u , and 
implies that there are many vacua that differ by their B, L content, separated by a potential 
barrier of electroweak height. At temperatures T «C 100 GeV transition between different vacua 
are negligible because suppressed by a quantum-tunneling factor e -2 ""/" 2 . If 2Vg en = 1 this would 
imply proton decay with an unobservably slow rate; since there are 3 generations and all of them 
must be involved, proton decay is kinematically forbidden. This suppression disappears at high 
temperature, T J> 100 GeV, and the space-time density of B, L-violating 'sphaleron' interactions 
is 7 ~ arf^T) faster than the expansion rate of the universe up to temperatures of about T ~ 
10 12 GeV [3]. 1 

3. SM baryogenesis is not possible due to the lack of out-of equilibrium conditions. The electroweak 
phase transition was regarded as a potential out-of equilibrium stage, but experiments now 
demand a higgs mass nih ^ 115 GeV, and SM computations of the Higgs thermal potential show 
that, for nih ^ 70 GeV, the higgs vev shifts smoothly from (h) = to (h) = v when the universe 
cools down below T ~ nih [3]- 

2. In any case, the amount of CP violation provided by the CKM phase would have been too small 
for generating the observed baryon asymmetry, because it is suppressed by small quark masses. 
Indeed CP-violation would be absent if the light quarks were massless. 

Many extensions of the SM could generate the observed n#. 'Baryogenesis at the electroweak 
phase transition' needs new particles coupled to the higgs in order to obtain a out-of-equilibrium 
phase transition and to provide extra sources of CP violation. This already disfavored possibility 
will be tested at future accelerators. 'Baryogenesis from decays of GUT particles' seems to conflict 
with non-observation of magnetic monopoles, at least in simplest models. Furthermore minimal GUT 
model do not violate B — L, so that sphaleron processes would later wash out the eventually generated 
baryon asymmetry. 

The existence of sphalerons suggests baryogenesis through leptogenesis: lepton number might be 
violated by some non SM physics, giving rise to a lepton asymmetry, which is converted into the 
observed baryon asymmetry by sphalerons. 

This scenario can be realized in many different ways [3]. Majorana neutrino masses violate lepton 
number and presumably CP, but do not provide enough out-of-equilibrium processes. The minimal 
successful implementation [4] needs just the minimal amount of new physics which can give the ob- 
served small neutrino masses via the see-saw mechanism [2]: heavy right-handed neutrinos N with 
masses M. 'Baryogenesis via thermal leptogenesis' [4] proceeds at T ~ M, when out-of-equilibrium 
(condition 3) CP-violating (condition 2) decays of heavy right-handed neutrinos generate a lepton 
asymmetry, converted in baryon asymmetry by SM sphalerons (condition 1). 

A real understanding of these issues needs advanced quantum field theory. This kind of theoretical studies lead to 
one observed consequence: the rj' mass, that is related to some QCD analogous of the SU(2)l effects we are considering. 
Therefore there should be no doubt that B, L are violated, and this is almost all what one needs to know to understand 
leptogenesis quantitatively. 
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Figure 2: CP-violating Ni decay. 

3 Thermal leptogenesis: the basic physics 

We now discuss the basic physics, obtaining estimates for the main results. The SM is extended by 
adding the heavy right-handed neutrinos suggested by see-saw models. To get the essential points, we 
consider a simplified model with one lepton doublet L and two right-handed neutrinos, that we name 
'Ni and l N 2j3 ', with Ni lighter than N 2j3 . The relevant Lagrangian is 

M 

X = JSfcM + iVi^iVi + AiiViiJL + -±Nl + (6) 
+N 2 , 3 i$N 2 , 3 + X 2 , 3 N 2<3 HL + -|^JVf, 3 + h.c. 

By redefining the phases of the N%, N 23 , L fields one can set Mi, M 23 , Ai real leaving an ineliminabile 
CP-violating phase in A2,3- 

3.1 CP-asymmetry 

The tree-level decay width of Ni is Ti = X 2 Mi/8tt. The interference between tree and loop diagrams 
shown in fig. 2 renders Ni decays CP-asymmetric: 



T(Ni — > LH) - T(Ni — > LH) 1 Mi - 

£l = — ; ; ; ^-=r- ~ lm A 9 o 

T(Ni -» LH) + r(JVi -» LH) 4vr M 2 , 3 2,3 

In fact 

T(Ni -> LH) oc |Ai + AX*i\l 3 \ 2 , T{Ni LH) oc |A* + ^4AiA|* 3 | 2 

where A is the complex CP-conserving loop factor. In the limit M 2j3 3> Mi the sum of the two 
one loop diagrams reduces to an insertion of the (LH) 2 neutrino mass operator mediated by ^2,3: 
therefore A is suppressed by one power of M 2)3 . The intermediate states in the loop diagrams in fig. 2 
can be on shell; therefore the Cutkosky rule guarantees that A has an imaginary part. Inserting the 
numerical factor valid in the limit M2.3 3> Mi we can rewrite the CP-asymmetry as 

3 Milmm 2 ,3 _ 6 Imm 2 , 3 Mi 

£l ~ 16vr v 2 0.05 eV 10 10 GeV [ ) 

where rh 2:3 = X 2 3 v 2 /M 2i3 is the contribution to the light neutrino mass mediated by the heavy ^2,3. 

The operator argument implies that eq. (7) holds in any model where particles much heavier than 
Mi mediate a neutrino mass operator with coefficient 1722,3. I n the past it was debated about if only 
the 'vertex' diagram in fig. 2 or also the 'self-energy' diagram should be included when computing the 
CP asymmetry: the operator argument makes clear that both diagrams contribute, since in the limit 
M2,3 3> Mi the two diagrams reduce to the same insertion of the (LH) 2 operator. 
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y NS = 7eq (LH -* LH) 



y m = 7eq (LL h> /f h; 



Figure 3: Wash-out LH <-» LLf and LL <-> LfLf AL = 2 scatterings. 



The final amount of baryon asymmetry can be written as 



n B 



9SM 



(8) 



where <7sm = 118 is the number of spin-degrees of freedom of SM particles (present in the denomi- 
nator of eq. (8) because only N\ among the many other particles in the thermal bath generates the 
asymmetry) and r/ is an efficiency factor that depends on how much out-of-equilibrium Aq-decays are. 



3.2 Efficiency 

We now discuss this issue. If N% — > LH decays are slow enough, the N\ abundancy does not decrease 
according to the Boltzmann equilibrium statistics un 1 oc e~ Ml ^ T demanded by thermal equilibrium, 
so that late out-of-equilibrium N\ decays generate a lepton asymmetry. Slow enough decay means N\ 
lifetime longer than the inverse expansion rate. At T ~ M\ one has 



R = 



r, 



H (M\ 



mi 

m* 



where 



m 



3M, 



2.3 W' 6 eV 



pi 



is fixed by cosmology. Ali the dependence on the mass and Yukawa couplings of Aq is incorporated in 
rhi = Afu 2 /Mi, the contribution to the light neutrino mass mediated by Aq. Unfortunately rh\ and 
7712,3 are only related to the observed atmospheric and solar mass splittings in a model-dependent way. 
Unless neutrinos are almost degenerate (and unless there are cancellations) rhi and 7712,3 are smaller 
than m atm ~ 0.05 eV. 

I£ R 1 (i. e. Ni decays strongly out-of-equilibrium) then rj = 1. 

If instead R^$> 1 the lepton asymmetry is only mildly suppressed as rj ~ l/R. The reason is that 
Ni inverse-decays, which tend to maintain thermal equilibrium by regenerating decayed Ni, have rates 
suppressed by a Boltzmann factor at T < M±: R(T < Mi) ~ R ■ e~ Ml ^ T . The Ni quanta that decay 
when R(T) < 1, i. e. at T < Mi/lnR, give rise to unwashed leptonic asymmetry. At this stage the 
Ni abundancy is suppressed by the Boltzmann factor e~ Ml ^ T = l/R. In conclusion, the suppression 
factor is approximately given by 



r] ~ min(l,m*/mi) (if Ni initially have thermal abundancy). 



(9) 



Furthermore, we have to take into account that virtual exchange of N\ t 2 t 3 gives rise to AL = 2 
scatterings (see fig. 3) that wash-out the lepton asymmetry. Their thermally-averaged interaction rates 
are relevant only at Mi > 10 14 GeV, when iVi 2,3 have large 0(1) Yukawa couplings. When relevant, 
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these scatterings give a strong exponential suppression of the baryon symmetry, because their rates 
are not suppressed at T<M\ by a Boltzmann factor (no massive N\ needs to be produced). 

So far we assumed that right-handed neutrinos have thermal initial abundancy. Let us discuss 
how the result depends on this assumption. If rh\ 3> m* (in particular if rh\ = m a t m or m sun ) the 
efficiency does not depend on the assumed initial conditions, because decays and inverse-decays bring 
the Nx abundancy close to thermal equilibrium. For fh\ < m* the result depends on the unknown 
initial condition. 

• If N\ have negligible initial abundancy at T > M\ and are generated only by the processes 
previously discussed, their abundancy at T ~ M\ is suppressed by rhi/rh*. Therefore the 
efficiency factor is approximatively given by 

r] ~ min(mi/m*, m* /fhi) (if N-y initially have zero abundancy). (10) 



• Finally, in the opposite limit where N\ initially dominate the energy density of the universe, 
the suppression factor 1/(/sm in eq. (8) no longer applies, and the efficiency factor can reach 
V ~ 5SM- 

r/ ~ min(<7sM5 rh*/rhi) (if iVi initially have dominant abundancy). (11) 



These estimates agree with the results of a detailed numerical computation, shown in fig. 4a. Notice 
that if fh\ 3> m* the N\ abundancy gets close enough to thermal equilibrium, such that the lepton 
asymmetry generated by N\ decays does not depend on the initial iVi abundancy. Furthermore, N\ 
decays and inverse-decays typically wash-out a possible pre-existing lepton asymmetry. 

In the next section we describe how a precise computation can be done. This is not necessary for 
understanding the final discussion of section 5. 



4 Thermal leptogenesis: precise computation 

As previously discussed many important computations in cosmology are done using Boltzmann equa- 
tions. So it is useful to have this tool. 



4.1 Boltzmann equations 

In absence of interactions the number of particles in a comoving volume V remains constant. Boltz- 
mann equations allow to follow the effect of different interactions. Let us study e. g. how 1 <-> 2 + 3 
decay and inverse-decay processes affect the number n\ of 'l' particles (in the case of leptogenesis we 
have N\ <-> LH): 



d 
di 



(niV) = V J dPi J dp2 J dp 3 (2TT) i 5\ Pl -p 2 - P3 ) x (12) 
x[-|^ 23 | 2 /i(l ± / 2 )(1 ± /s) + 1^23-1 1 2 (1 ± /l)/ 2 / 3 ] 



where dpi = d 3 pj/2£'j(27r) 3 is the relativistic phase space, |Ai_>23| 2 and 1^423^1 1 2 are the squared 
transition amplitude summed over initial and final spins, and f i are the energy (and eventually spin, 
flavour, color,...) distributions of the various particles. To start we assume that CP violation can be 
neglected, such that the direct and the inverse process have a common amplitude A. 

In line of principle we should study the evolution of all /'s in order to obtain the total densities 
n = ^2 f f d 3 p/(2ir) 3 . In practice elastic scatterings (i. e. interactions that do not change the number of 
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Figure 4: Fig. 4a: efficiency rj as function of rh\ for Mi = 10 10 GeV and for different assumptions 
about the initial N\ abundancy; for fh\ larger than a few meV the efficiency is univocally predicted to 
be 77(7711) ~ 0.42(meV/mi) 1 ' 15 . Fig. 46: the regions in the (777,1, Mi) plane inside the curves can lead 
to successful leptogenesis. 

particles) are typically fast enough that they maintain kinetic equilibrium, so that the full Boltzmann 
equations for / are solved by f (p) = f C q(p)n j 'n cq where / cq = [e E ' T ± are the Bose-Einstein and 
Fermi-Dirac distributions. Each particle species is simply characterized by its total abundancy n, that 
can be varied only by inelastic processes. 

The factors l±/j in eq. (12) take into account Pauli-Blocking (for fermions) and stimulated emission 
(for bosons). Since the average energy is {E} ~ 3T within 10% accuracy one can approximate with 
the Boltzmann distribution / cq ~ e~ E l T and set 1 ± / ~ 1. This is a significant simplification. 

When inelastic processes are sumciently fast to maintain also chemical equilibrium, the total num- 
ber n eq of particles with mass M at temperature T are 

d 3 P / cq _ gM 2 T M f gT 3 /7T 2 T > M 



Heq ~ 9 J (2TThf ~ 2vr2 K2( T ] ~ \ g(MT/27tf/ 2 e- M / T T « M 

where g is the number of spin, gauge, etc degrees of freedom. A right handed neutrino has g^ = 2, 
a photon has <? 7 = 2, the 8 gluons have gc a = 16, and all SM particles have #sm — 118. The factor 
h = h/2ir has been explicitly shown to clarify the physical origin of the 2ir in the denominator. 
The Boltzmann equation for 77-1 simplifies to 

V~d~t^ niV ^ = J d ^ J d ^ 2 J d P3( 27r ) 4 ^iPi-P2-P3) x (13) 
AA?[-^e-^ T + ^e-^ T e-^] 

One can recognize that the integrals over final-state momenta reconstruct the decay rate Vi, and that 
the integral over d?p\/ E\ averages it over the thermal distribution of initial state particles; the factor 
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1/Ei corresponds to Lorentz dilatation of their life-time. Therefore the final result is 



= < r ^ 



n 



cq 



n 2 n 3 

eq eq 

n n 



v l ,b 2 "3 

where (ri) is the thermal average of the Lorentz-dilatated decay width 



ri(£l) = ^~J d P2 dp 3 (27T)*5 4 ( Pl -P2- P3 )\A\ 2 . 



(14) 



(15) 



Analogous results holds for scattering processes. 

If (ri) 3> H the term in square brackets in eq. (14) is forced to vanish. This just means that 
interactions much faster than the expansion rate force chemical equilibrium, giving n = n eq . 

In the case of leptogenesis 2,3 = L, H have fast gauge interactions. Therefore we do not have to 
write and solve Boltzmann equations for L, H, because we already know their solution: L, H are kept 
in equilibrium. We only need to insert this result in the Boltzmann equation for N\, that simplifies to 



h\ + 3-ffni = (ri)(ni - n^) 



(16) 



having used V /V = 3H = —s/s. 

In computer codes one prefers to avoid very big or very small numbers: it is convenient to reabsorb 
the 3H term (that accounts for the dilution due to the overall expansion of the universe) by normalizing 
the number density n to the entropy density s. Therefore we study the evolution of Y = n/s, as 
function of z = m^/T in place of time t (H dt = d\nR = dlnz since during adiabatic expansion sV 
is constant, i. e. V oc l/T 3 ). 

Using Y (z) as variables, the general form of Boltzmann equations is 



dz 



34- 



J3_Y4_ 
z 3 J 4 



(17) 



where the sum runs over all processes that vary the number of 'l' particles by Ai units (e. g. Ai = — 1 
for 1 — > 23 decay Ai = —2 for 11 — > 23 scatterings, etc.) and 7 eq is the specetime density (i. e. the 
number per unit volume and unit time) in thermal equilibrium of the various processes. 

Neglecting CP-violating effects, direct and inverse processes have the same reaction densities. For 
a scattering and for its inverse process one gets the previous result: 



7eq (l - 23) = / dpx fl 



r><i 



93 (27r)^( Pl -p 2 - p 3 )\A\ 2 = 7eq (23^ 1) 
The thermal average of the decay rate can be analytically computed in terms of Bessel functions: 

eqKl(z) 



7eq(23' 



D 



1 K 2 (z) 



r(l -> 23' 



7eq (l 23 ■ ■ ■) 

For a 2-body scattering process 

7cq (12 ^ 34) = y dpi dp 2 /r/ 2 Cq / dp 3 dp A (2tt) 4 <5 4 (pi + p 2 - Ps - p A )\A\' 



(18) 



When there are n identical particles in the final or initial states one should divide by a n! symmetry 
factor. One can analytically do almost all integrals, and obtain 



7eq (12->34"0 



32tH 



ds s 3/2 \(l,Mf/s,M$/s)a(s) Kj 



T 



which is the thermal average of v ■ a, summed over initial and final state spins. 
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Figure 5: Evolution ofY^^Mi/T) for different values of rh\ starting from different initial conditions. 
The continuous line is the abundancy in thermal equilibrium. 



4.2 Boltzmann equations for leptogenesis 

We now specialize to leptogenesis. The main process is N\ <-> H L, H L decay and inverse decay. These 
processes are enough to generate a N± population that follows the Boltzmann distribution 2 , so that 
we can write the Boltzmann equation for the total N\ abundancy. We denote with 7^ its equilibrium 
density rate, computed inser ting r(iVi) = AfMi/8-Tr in eq. (18). The Boltzmann equation for the N\ 
abundancy is 

dY Nl 



sHz- 



dz 



D. 



(19) 



Fig. 5 shows howe YXr 1 evolves for different values of rh\. As expected if rh\ S> 10~ 3 eV one gets a 
result close to thermal equilibrium independently of the initial condition. 

In order to get the Boltzmann equation for the lepton asymmetry one needs to take into account 
the small CP-violating terms. Let us start by including only the AL = 1 CP-violating Ni — > H L, H L 
decays. We write the decay rates in terms of the CP-conserving total decay rate 7d and of the 
CP-asymmetry e\ <C 1: 



7cq (iVi - LH) = 7eq (LH - Ni) = (1 + £i)7d/2, 
7eq (iVi -> LH) = leq (LH Nx) = (1 - ei) 7D /2. 



In this approximation the number of leptons L and anti-leptons L evolve as 



zHsY' L 



1D 

2 



— (1 



r; 



(20) 



(21) 



2 Intuitively, one would probably guess a different, incorrect, result: that H L — > JVi inverse decays generate JVi with the 
average energy of two particles, rather than with the average energy of one particle, since Em x — El + Eh, and that two 
particles have more energy than one. However, using e~ Ei//T ■ e _Eff//T = e~ BjY i'' T one verifies that (En-l) = (El + Eh): 
this is why Boltzmann found that the thermal distribution is exponential. 
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zHsYl = ^ 

L o 



^(i_ ei )_^(i + ei ; 



(22) 



Here Y^, Y^ q and Y^ q are equilibrium densities each with 2 degrees of freedom. Ignoring O(ef) terms, 
the lepton asymmetry C = Yl — Yi evolves as 

N L1 L 



sHzC' = £i7d(S + 1) - J^lD- (23) 



The second term describes how 7/5 tends to restore thermal equilibrium, washing out the lepton 
asymmetry. The first term makes no sense: it would generate a lepton asymmetry even in thermal 
equilibrium, Yn 1 = Y^, violating Sakharov conditions. An acceptable Boltzmann equation would 
contain Y^/Y^ — 1, but we made no sign error. Indeed, taking into account only decays and inverse 
decays, an asymmetry is really generated even in thermal equilibrium, since CPT invariance implies 
that if iVi decays preferentially produce L, than inverse decays preferentially destroy L i. e. they have 
the same net effect. 



A subtlety: avoiding overcounting 

To obtain correct Boltzmann equations one must include all processes that contribute at the chosen 
order in the couplings. The CP-asymmetry is generated at O (A 4 ): at this order we must include also 
the AL = 2 scatterings of fig. 3; we name their rates as 

INs = J eq (LH LH), lNt = 7cq (LL H H), (24) 

and 7al=2 = 'ZiiNs+iNt)- At first sight it is enough to include these scatterings at tree level, obtaining 
CP-conserving reaction densities 7 = 0{\ ) that cannot correct our non-sensical CP-violating term. 
This is basically right, although the true argument is more subtle. 

Indeed LH <-> LH can be mediated by on-shell N\ exchange (see fig. 3a): as usual in these 
situations (e. g. the Z-peak) resonant enhancement gives cr pca k oc A° in an energy range AE oc A 2 , 
so that 7tv s oc cr pca k • AE oc A 2 . (We will soon obtain the exact result). Nevertheless one can proof 
that the reaction density remains CP-conserving up to one-loop order: unitarity demands ^ ■ \M(i — > 

i)l a = £; W-0l a .*> 

o(LH -> LH) + o(LH -► LH) = a{LH -» LH) + a(LH -> LH) 

(at higher order states with more particles allow a negligible CP asymmetry). 

The key subtlety is that the LH <-> LH scattering rate mediated by s-channel exchange of N\ 
shown in fig. 3a, must be computed by subtracting the CP-violating contribution due to on-shell N\ 
exchange, because in the Boltzmann equations this effect is already taken into account by successive 
decays, LH <-> N\ <-» LH. Since the on-shell contribution is 

^- shcU (LH - LH) = y cq (LH - N 1 )BR(N 1 - LH), 

where BR(Aq -> LH) = (1 - £i)/2, we obtain 

y™ h (LH^LH) = lNs -(l- £l f lD /4, (25) 
^(LH^LH) = lNs -(l +£l ) 2 lD /4, (26) 

Including subtracted scatterings at leading order in £\ gives the final Boltzmann equation [4] 



zHsC' = 7d 



2 7AL=2yeq 5 ( 27 ) 
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where 7al =2 = 7al=2 - 7d/ 4 - 

Equivalently, one can more simply not include the decay contribution 70 to washout of £ because 
it is already accounted by resonant AL = 2 scatterings. Then one gets 

zHsC = 7D ei ( ^ - 1 j - 2 7AL=2 A (28) 



which is equivalent to (27). Using eq. (28) it is not necessary to compute subtracted rates. (Subtraction 
is performed incorrectly in works before 2003). 

Including sphalerons and Yukawas 

Finally, we have to include sphaleronic scatterings, that transmit the asymmetry from left-handed 
leptons to left-handed quarks, generating a baryon asymmetry. Similarly we have to include SM 
Yukawa couplings, that transmit the asymmetry to right-handed quarks and leptons. 

In theory one should enlarge Boltzmann equations adding all these processes. In practice, depend- 
ing on the value of M\, during leptogenesis at T ~ M\ these process often give reaction densities 
which are either negligibly slower or much faster than the expansion rate: in the first case they can 
be simply neglected, in the second case they simply enforce thermal equilibrium. One can proceed 
by converting the Boltzmann equation for C into a Boltzmann equation for B — C: since B — C is 
not affected by these redistributor processes we only need to find how processes in thermal equilib- 
rium relate B — C to C. Sphalerons and the Xt,b,c,r Yukawas are fast at T <^ 10 n ~ 12 GeV. At larger 
temperatures all redistribution processes are negligibly slow and one trivially has B — C = —C. At 
intermediate temperatures one has to care about flavor issues, discussed later. 

It is interesting to explicitly compute redistribution factors at T ~ TeV when all redistributor 
processes are fast. Each particle P = {L, E, Q, U, D, H} carries an asymmetry Ap. Interactions 
equilibrate 'chemical potentials' fip = Ap/gp as 

ELH Yukawa : = (j,e + I^l + /J- H 

DQH Yukawa : = fip, + (j,q + fijj 

U QH Yukawa : = fijj + /ig — fin 

QQQL sphalerons : = 3[iq + /j-l 

No electric charge : = N gcn (fj, Q - 2/ij/ + /iD — I^L + He) - 2N mggs n H 

Solving the system of 5 equations and 6 unknowns, one can express all asymmetries in terms of one 
of them, conveniently chosen to be B — C: 

28 

B = N gen (2fi Q - fiu - fi D ) = —(B - C), 

^1 

C = B-(B-C) = --(B-C). 
The efficiency r\ is precisely defined such that the final baryon asymmetry is 

etrj 



- = B = --e V Y N «(T » Ml ) ,e. - 



today 103- 



(29) 



in agreement with the estimate (8). 

Various extra processes give corrections of relative order g 2 /n 2 , A 2 /V 2 ~ few %. Some of these 
corrections have been already computed: scattering involving gauge bosons and/or top quarks. Others 
have not yet been included: three body iVi-decays, one-loop correction to the Ni — > LH decay and 
its CP-asymmetry. Thermal corrections have not been fully included. The fact that 7^ is the only 
really relevant rate makes a full inclusion of these subleading corrections feasible. 
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Including flavor 

So far we studied the dynamics of leptogenesis in 'one-flavor' approximation, eq. (6). In the literature, 
flavor was fully included only after that neutrino data showed that flavor mixing among leptons can be 
large. The one-flavor approximation remains useful because including flavor adds so many unknown 
parameters that a precise discussion is impractical: e. g. we do not know which combination of L e , L„, 
L T is the lepton doublet L coupled to N\ in the see-saw Lagangian, eq. (6). To include flavor, the 
Bolztmann equation for Yc must be generalized to an evolution equation for the 3x3 matrbc density 
p of lepton asymmetries in each flavor. However, it simplifies to qualitatively different behaviors in 
different ranges of M\\ 

• If Mi J> 10 11 GeV all SM lepton Yukawa couplings induce scattering rates much slower than the 
expansion rate H at T ~ M\: quantum coherence among different flavors survives undamped 
and the main new effect is that lepton asymmetries generated by N2 , N% decays can be washed- 
out by processes involving N\ only along the combination of flavors to which N\ couples. So, 
one must sum the contributions of all right-handed neutrinos produced after reheating. 

• If M\ < 10 9 GeV, the and A r Yukawa couplings induce scattering rates faster than H at 
T ~ Mi and damp quantum coherence in p: the matrhc equation for p reduces to 3 Boltzmann 
equations for the asymmetries in the £ = {e, p, r} flavors. Neglecting a mild flavor mhdng 
(induced by sphaleronic scatterings) these equations have the following form: eq. (19) for Yj^ 
remains unchanged, and eq. (28) for C splits into three equations for C e , £ M , C T with the CP- 
violating term and the wash-out term restricted to each flavor, as intuitively expected. Therefore, 
in each flavor l one has a different CP asymmetry E\i and efficiency T]t{jh\^rh\i) that now depends 
on two parameters: the usual flavor independent rh\ that tells the total N\ decay rate, and the 
flavor-dependent rhi e u,r = \^u\ v2 /M± that parameterize wash-out. A good approximation is [4]: 

~ r](m u ). (30) 

where 77 is the one-flavor efficiency, plotted in fig. 4a. The 
side figure shows the exact numerical result of r]£ as function 
of rhu for different values of fh\/rh\i = 1 (continuous line), 
2 (red dotted line), 4 (blue dashed line), confirming that the 
result negligibly depends on mi, especially if it is close to the 
values suggested by solar or atmospheric oscillations. There- 
fore eq. (8) gets replaced by 

ub/ti^ & y~)eu ■ r}{m u ) / gsu- (31) 
1 

For large mhdng angles one typically has rhi/rhu ~ 0(2): 
the above approximation tells that taking flavor into account 
can enhance the efficiency by an O (2) factor. 

• Something intermediate happens if 10 9 GeV < Mi < 10 11 GeV: quantum coherence stays un- 
damped only among p and e, such that one must take into account the asymmetry possibly 
generated by N2, N3 along e, p. 




w 5 w 4 i(r 3 10 2 10-' 1 



m 1 f of wash - out in eV 



14 



5 Testing leptogenesis? 



Unfortunately speculating that neutrino masses and the baryon asymmetry are produced by the see- 
saw mechanism and by thermal leptogenesis is much easier than testing them. 

A direct test seems impossible, because right-handed neutrinos are either too heavy or too weakly 
coupled to be produced in accelerators. 

What about indirect tests? We trust BBN because it predicts the primordial abundances of several 
nuclei in terms of known particle physics. Leptogenesis explains a single number, ns/n 7 , in terms 
of speculative physics at high energies: the see-saw model has 18 free parameters. Neutrino masses 
only allow to measure 9 combinations of these parameters, and thereby provide a too weak link. 
The situation might improve if future experiments will confirm certain supersymmetric models, in 
which quantum corrections imprint neutrino Yukawa couplings XijNiLjH in slepton masses, inducing 
lepton flavour violating (LFV) processes such as \i — > ej, r — ► //7, r — ► e7 with possibly detectable 
rates. Measuring them, in absence of other sources of LFV, would roughly allow us to measure 
the 3 off-diagonal entries of A^ • A. Detectable LFV rates are obtained if A> KH 1 " 2 ). In any case, 
reconstructing all see-saw parameters in this way is unrealistic. Maybe future experiments will discover 
supersymmetry, LFV in charged leptons, and will confirm that neutrino masses violate lepton number 
and CP, and we will be able to convincingly argue that this can be considered as circumstantial 
evidence for see-saw and thermal leptogenesis. Archeology is not an exact science. 

Another possibility is that we might find a correctly predictive model of flavour. Presently three 
approaches give some predictions: symmetries, numerology, zerology. Symmetries can be used to 
enforce relations like #23 = vr/4, tan 2 #12 = 1/2, #13 = where Oij are the neutrino mixing angles. 
Numerology can suggest relations like 9\2 + Oc = tt/4. Zerology consists in assuming that flavour 
matrices have many negligibly small entries; for example one can write see-saw textures with only 
one CP-violating phase. This scheme does not allow to predict the sign of CP-violation in neutrino 
oscillations in terms of the sign of the baryon asymmetry, because the sign of the baryon asymmetry 
also depends on which right-handed neutrino is the lightest one. 



5.1 Constraints from leptogenesis 

We here discuss testable constraints. Although this topic is tortuous, we avoid over-simplifications, at 
the price of obtaining a tortuous section. 

As discussed in section 3, assuming that N\ is lighter enough than other sources of neutrino masses 
that their effects can be fully encoded in the neutrino mass operator, the CP-asymmetry s\ is directly 
connected with it. Under the above hypothesis, one can derive constraints from leptogenesis. To do so, 
we need to generalize eq. (7) taking flavour into account. Denoting flavour matrices with boldface, we 
define m; to be the contribution to the neutrino mass matrbc mediated by the right-handed neutrino 
Ni, so that m u = rh\ + rfi2 + rh^. Then 

N = 3* ts^!i^ < 3 jft 

l07T V z mi l07T V A 

where m„ 3 (m ui ) denotes the mass of the heaviest (lightest) neutrino. Rather than rigorously proofing 
the constraint (32) [5], let us understand its origin and limitations in a simpler way. 

1) Let us start considering the case of hierarchical neutrinos: m Ul «C m V3 : the constraint is ob- 
tained by substituting rh,2 + rh^ = m v — rhi, and holds whatever new physics produces m2 3- 
Since |ei| oc Mi, one can derive a lower bound on the mass Mi of the lightest right-handed 
neutrino by combining eq. (32) with a precise computation of thermal leptogenesis in one-flavor 
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approximation and with the measured baryon asymmetry and neutrino masses: 

( 2.4 x 10 9 GeV if Ni has zero 
Mi > < 4.9 x 10 8 GeV if Ni has thermal initial abundancy (33) 
{ 1.7 x 10 7 GeV if Ni has dominant 

and assuming Ml <C M23. Including flavor as discussed previously relaxes the first constraint 
by a factor O (2). 

2) The factor m Vz — m Ul is specific to the see-saw model with 3 right-handed neutrinos. To under- 
stand it, let us notice that 3 right-handed neutrinos can produce the limiting case of degenerate 
neutrinos m Ul = m U2 = m V3 = m v only in the following way: each Ni gives mass m„ to one neu- 
trino mass eigenstate. Since they are orthogonal in flavour space, the CP-asymmetry of eq. (32) 
vanishes due to flavour orthogonality: this is the origin of the m U3 — m Ul suppression factor. 

The bound (32) implies an upper bound on the mass of quasi-degenerate neutrinos: m„ < 0.2 eV [5] 
Indeed for large m v both the efficiency and the maximal CP-asymmetry become smaller, because 
heavy neutrinos must be quasi-degenerate, and rh\ > m ui , m„ 3 — m Vl ~ Am 2 tm /2m l ,. Further- 
more, the bound can be improved up to m Us < 0.15 eV [5] (3<r confidence level) by computing 
the upper bound on |ei | for given fh\ and maximizing n b with respect to fh\ taking into account 
how the efficiency of thermal leptogenesis decreases for large rh\. The leptogenesis constraint is 
very close to observed neutrino masses, and stronger than experimental bounds. 

However, this leptogenesis constraint holds under the dubious assumption that hierarchical right- 
handed neutrinos produce quasi-degenerate neutrinos, while good taste suggests that quasi-degenerate 
neutrinos are more natur ally produced by quasi-degenerate right-handed neutrinos. In general the 
constraint (32) evaporates if the particles that mediate /77,23 are so light that their effects cannot 
be encoded in JTI23: the CP-asymmetries becomes sensitive to the detailed structure of the neutrino 
mass model. Suppression due to flavor-orthogonality was the most delicate consequence of our initial 
assumption, and is the first result that disappears when they are relaxed, allowing M2,3 to be not 
much heavier than Mi. Correspondingly, the constraint on quasi-degenerate neutrino masses, that 
heavily relies on the factor m U3 — m Ul , becomes weaker than experimental constraints. 

Furthermore, in the extreme situation where right-handed neutrinos are very degenerate, M2 — 
Mi ^Ti^ a qualitatively new effect appears: CP violation in Ni «-> N2 mixing. This phenomenon is 
fully analogous to K° «-> K° mixing, and for M2 — Mi ~ Fi 2 it allows a maximal CP-asymmetry, 
|ei| ~ 1. This means that with a tiny M2 — Mi one can have successful thermal leptogenesis even at 
the weak scale. 

The constraint (33) is more robust, but we still have to clarify what the assumption M2,3 3> Mi 
means in practice. Surely one needs M2 — Mi 3> ri 2 such that only CP-violation in N\ decay is 
relevant. The issue is: is M2$/M\ ~ 10 (a hierarchy stronger than the one present in left-handed 
neutrinos) hierarchical enough to guarantee that the constraint holds, up to 1/10 2 = % corrections? 
The answer is no: an operator analysis allows to understand how the constraint can be completely 
relaxed. The physics that above Mi produces the dimension-5 neutrino mass operator (LH) 2 /2 can 
also produce a related dimension-7 operator T = (LH)d 2 (LH)/2, that contributes to E\. Since 
neutrino masses do not constrain T, it can be large enough to over-compensate the suppression due 
to its higher dimension. We make the argument more explicit, by considering the concrete case of 
see-saw models, where above Mi there are two other right-handed neutrinos with masses M2 and 
M3. Including 'dimension-7' terms suppressed by Mf/M| 3 and dropping inessential 0(1) and flavour 
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factors, the CP-asymmetry becomes: 



3 Mi 

167T f 2 



Im 



~ /, M f\ ~ , n M i- 
m2(1 + Mf )+m3(1 + M|' 



(34) 



At leading order e\ depends only on rfi2 + A3, that, as previously discussed, cannot be large and 
complex. At higher order e\ depends separately on m,2 and 77J3, that can be large and complex 
provided that their sum stays small. One can build models where this naturally happens, obtaining 
an ei orders of magnitudes above the DI bound. The enhancement is limited only by perturbativity, 
^2,3 ^47r. In supersymmetric models, large Yukawa couplings lead to the testable effects previously 
discussed. 



5.2 Leptogenesis and supersymmetry 

Adding supersymmetry affects some 0(1) factors. Eq. (29) remains almost unchanged, because adding 
spartners roughly doubles both the number of particles that produce the baryon asymmetry and 
the number of particles that share it. Ignoring small supersymmetry-breaking terms, right-handed 
neutrinos and sneutrinos have equal masses, equal decay rates and equal CP-asymmetries. Both 
and Ei become 2 times larger, because there are new decay channels. As a consequence of more CP- 
asymmetry compensated by more wash-out, the constraints on right-handed and left-handed neutrino 
masses discussed in the non-supersymmetric case remain essentially unchanged. 

The leptogenesis constraint (33) on Mi acquires a new important impact: in many supersymmetric 
models the maximal temperature at which the Big-Bang started (or, more precisely, the 'reheating tem- 
perature' Trh) must be less than about Trh ^ 10 GeV, in potential conflict with the constraint (33). 

Let us discuss the origin of the supersymmetric constraint on Trh- Gravitinos are the supersym- 
metric partner of the graviton: they have a mass presumably not much heavier than other sparticles, 
m,Q < TeV, and gravitational couplings to SM particles. Therefore gravitinos decay slowly after BBN 
(life-time r ~ Mpjm^, ~ sec (lOOTeV/m^) 3 ): their decay products damage the nuclei generated 
by BBN. The resulting bound on the gravitino abundancy depends on unknown gravitino branching 
ratios. 3 The gravitino interaction rate is Jq{T) ~ T^/Mpj, which means that gravitinos have been 
generated around the reheating temperature Irh, with abundancy n^/n^ ~ j^/Hn^ ~ T^n/Mpi. 
Therefore gravitinos suggest an upper bound on T. 

This is why in the previous section we carefully discussed specific scenarios that allow low-scale 
thermal leptogenesis. Supersymmetry suggests a new scenario named 'soft leptogenesis': complex soft 
terms in the see-saw sector give new contributions to the CP-asymmetry, z\ ~ a2 m -|usY/-^i ' which 
can be significant if M\ is not much larger than the SUSY-breaking scale msusY (presumed to be 
below 1 TeV). At larger M\ 'soft leptogenesis' can still be relevant, but only in a fine-tuned range of 
parameters. 



5.3 Leptogenesis in alternative neutrino-mass models 

Generic neutrino masses can be mediated by tree-level exchange of three different kinds of new par- 
ticles [2]: I) at least three fermion singlets; II) at least three fermion SU(2)i triplets; III) one scalar 
SU(2)^ triplet (or of combinations of the above possibilities). Fig. 6 shows the relevant Feynman 
diagrams. 

So far we studied case I). Can leptogenesis be used to distinguish between these possibilities? 
Leptogenesis can be produced in decays of P, the lightest particle that mediates neutrino masses. The 

3 Gravitinos might be stable if they are the lightest supersymmetric particle. But in this case dangerous effects are 
produced by gravitational decays of the next-to-lightest SUSY particle into gravitinos. 
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Figure 6: The Majorana neutrino mass operator (LH) 2 can be mediated by tree level exchange of: I) 
a fermion singlet; II) a fermion triplet; III) a scalar triplet. 

neutrino- mass contribution to its CP-asymmetry is given by expressions similar to eq. (32), with rh\ 
generalized to be the contribution to neutrino masses mediated by P, and rh,2 + rhs generalized to be 
the contribution of all heavier sources. 

It was expected that the efficiency r/ can be high enough only if P is a right-handed neutrino: as 
discussed in section 3.2 it easily decays out-of-equilibrium giving 



rj (fermion singlet) ~ min 



' r 



(35) 



where T is its decay rate, H is the expansion rate at T ~ M and 



{1 for thermal 

T/H for negligible initial abundancy. (36) 
<7SM for dominant 

A SU(2)i triplet (scalar or fermion) has gauge interactions that keep its abundancy close to thermal 
equilibrium so that the 3rd Sakharov condition cannot be fulfilled. This suppression is present, and 
a quantitative analysis is needed to see if/when it is strong enough. The Boltzmann equation for the 
triplet abundancy Y has an extra term 'ja that accounts for annihilations of two triplets into gauge 
bosons: 

dY Y Y 2 

sHz— = - 1D (— - 1) - 2 7A (-^ - 1). (37) 

1 e q 1 e q 

The term 7,4 ~ <? 4 T 4 is dominant only at T J> M, where M is the triplet mass; at lower temperatures it 
is strongly suppressed by a double Boltzmann factor (e -M / T ) 2 , because gauge scattering must produce 
2 triplets. The resulting efficiency can be approximated as 



//(fermion triplet) ~ min 



H M .r. 

''T'To^GeV 111 ^ 1 '^ 



(38) 



rj is univocally predicted, because at T 3> M gauge interactions thermalize the triplet abundancy. At 
T ~ Mi gauge scatterings partially annihilate triplets: in section 1.1 we learnt how to estimate how 
many particles survive to annihilations, and this is the origin of the factor M/10 12 GeV in rj. The last 
factor takes into account that annihilations are ineffective if triplets decay before annihilating. 

As illustrated in fig. 6 a scalar triplet separately decay to leptons, with width Tl(T* — > LL), and 
in Higgses, with width T h (T — > H H). Lepton number is effectively violated only when both processes 
are faster than the expansion of the universe, giving 



r/(scalar triplet) w min 



1 H M n r L+r H . 

' min(r L ,r^) ' 10 12 GeV maXl ' H ' 



(39) 
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This means that a quasi-maximal efficiency rj ~ 1 is obtained when T* — > LL is faster than gauge 
annihilations while H — > TT is slower than the expansion rate. In conclusion, leptogenesis from decays 
of a SU(2)x, triplet can be sufficiently efficient even if triplets are light enough to be tested at coming 
accelerators, M ~ TeV. 

5.4 Leptogenesis and dark matter in loop-mediated neutrino-mass models 

Neutrino masses mediated at one-loop level can be realized by exchanging various possible kinds of 
new particles; some of them are potential DM candidates. For example, one can introduce the usual 
right-handed neutrinos N and couple them as A N LH' , where H' is not the usual Higgs scalar doublet, 
but another scalar doublet coupled to the Higgs as X'(H^ H') 2 + h. c. The Lagrangian is restricted to 
couplings invariant under the Z2 symmetry N — > —N and H' — ► —H' , such that the lightest component 
of these particles is stable and is a DM candidate. 

Neutrino masses are suppressed by a one loop factor r ~ A' / (4vr) 2 with respect to the Standard 
see-saw. Leptogenesis depends on the coupling A and on M\ t 2,3 in the usual way: therefore, at fbced 
values of neutrino masses, the CP asymmetry (oc rfi2,3, see eq. (7)) is enhanced by l/r with respect to 
the Standard case, while the efficiency factor r/ (that depends on mj as in eq.s (9-11)) is suppressed 
by r if fh\ ~ m S un,atm- The constraints in eq. (33) on M\ are relaxed by the factor l/r. 

DM might be generated in the usual way discussed in section 1.1. Alternatively, a new interesting 
possibility is that both DM and leptogenesis might be generated by the same out-of-equilibrium process. 
In the example above, N% decays generate an asymmetry both in L and H': H' could be DM with 
abundancy dominated by its asymmetry (like protons) rather than by the usual freeze-out relic abun- 
dancy discussed in section 1.1. //the DM asymmetry is not washed-out, this scenario leads to a precise 
testable prediction for the DM mass: jtidm = cm p f2DM/^b ~ lOGeV, where c is a model-dependent 
O(l) factor, analogous to the 28/79 in eq. (29). 

However the DM asymmetry survives only until H gets a vev breaking the electroweak gauge 
symmetry: this happens at temperatures T < 1.2mh- Indeed DM is neutral under unbroken electro- 
magnetism and a Z2-like symmetry guarantees DM stability but does not protect the DM asymmetry. 
One can avoid wash-out by promoting Z2 to a global U(l): we now discuss why this class of models 
seems not viable. In the model used as concrete example, imposing a global U(l) means setting A' = 0, 
such that H' is not affected when H develops its vev. However, direct searches for DM told that DM 
must not only be neutral under electromagnetism, but also almost neutral under the Z: the DM Z 
coupling must be 10 -f- 100 times smaller than a typical electroweak coupling [6]. If the DM multiplet 
(either scalar of fermionic) lies in a complex representation of the electroweak gauge group (such that 
it can carry an asymmetry), reduced Z couplings need A'-like couplings. E. g., in the concrete H' 
example, A' splits the complex neutral component of H', S + iA, into a real scalar and pseudo-scalar, 
S and A, with different masses. The Z^ does not couple to the lightest DM mass eigenstate S or A but 
only to the off-diagonal combination A ■ d^S — S ■ d^A, without giving unseen DM signals if the mass 
difference Am = \ms — tua\ is bigger than about 20keV (the expected kinetic energy of DM around 
the earth) . One can show that the needed Am is so big that oscillations among S + i A and S — i A 
destroy the DM asymmetry. 4 To avoid this conclusion one needs either tbdm ~ SOrrih or models where 
DM is a complex neutral scalar singlet coupled to the Higgs: both possibilities look unattractive. 

Acknowledgments The original work presented in the last section was performed in collaboration 
with M. Raidal and V. Rychkov, that agreed on briefly summarizing here our no-go result. I thank 

4 This needs a non trivial analysis: Bolztmann equations for the density matrix of S and A show that DM-number- 
conserving gauge scatterings with rate r 3> H synchronize oscillations among different moment and reduce their frequency 
from Am to Am 2 /T. These non-trivial features help a lot, but not enough. 
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